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RESUMEN 

En este papel aplicamos y discutimos un metodo para la determination de la 
evolution del campo magnetico en los remanente de la supernova (SNRs) de la luminosi- 
dad de radio en la frecuencia v dada a la correlation del diametro (L„ — D). Asumimos 
que H se desarrolla como H oc D~ s , donde D estael diametro del remanente. El valor 
8 ~ 1.2 se obtiene bajo asuncion del equipartition de las ecuaciones para el calculo re- 
visado del equipartition (REC) y usar la muestra de los datos de la galaxia proxima M82 
del starburst. Intentamos investigar si o no los remanente de la supernova en M82 est 
en el estado del equipartition. Inconsistence entre el valor obtenido teoricamente para las 
condiciones y la extension adibatica (8 = 1.5) del equipartition y el valor empico obtenido 
adjunto, se puede explicar principalmente por la influencia de los efectos de la selection 
de la sensibilidad que tienden para aplanar la cuesta de las relaciones del L v — D para las 
muestras extragalactic. 

ABSTRACT 

In this paper we apply and discuss a method for the determination of the mag- 
netic field (H) evolution in supernova remnants (SNRs) from radio luminosity at given 
frequency v to diameter (L„ — D) correlation. We assumed that H evolves as H oc D~ s , 
where D is the diameter of the remnant. Value 8 « 1.2 is obtained under the equipartition 
assumption from the equations for revised equipartition calculation (REC) and by using 
the data sample from the nearby starburst galaxy M82. We try to investigate whether 
or not SNRs in M82 are in the equipartition state. This is done by comparison of our 
empirically obtained 8 with the theoretical value expected for the equipartition conditions. 
Inconsistence between the value obtained for equipartition conditions and adiabatic ex- 
pansion (S = 1.5) and the value empirically obtained herein, can be explained mainly by 
the influence of sensitivity selection effects which tend to flatten the slope of the L v — D 
relations for extragalactic samples. 

Key Words: ISM: MAGNETIC FIELDS — GALAXIES: INDIVIDUAL (M82) — 
METHODS: STATISTICAL — RADIATION MECHANISMS: NON- 
THERMAL — RADIO CONTINUUM: ISM — SUPERNOVA REM- 
NANTS 



1. INTRODUCTION 

Supernova remnants (SNRs) stand an important 
factor in the process of cosmic ray acceleration and 
matter circulation. Albeit very important, these pro- 
cesses are still not fully understood. Various theo- 
ries were suggested during the last few decades with 
a view to understanding the SNR properties. There 
is a general belief that the evolution of an SNR is 
strongly influenced by the properties of the local in- 
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terstellar medium (ISM) in which SNR evolves. As 
SNRs are the luminous synchrotron emitters in radio 
domain of the electromagnetic spectrum, the mag- 
netic field inside them and energetic spectrum of 
relativistic particles can be determined. Here, we 
will mainly focus on the magnetic field properties 
such as field strength and evolution. The most com- 
monly used empirical relation in studies of SNR evo- 
lution properties is radio surface brightness to di- 
ameter (£ — D) relation. This is because the only 
statistically reliable data samples of SNRs are found 
in radio domain. In order to study SNR evolution 
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issues from a slightly different perspective, in this pa- 
per we apply a method that transforms S — D into 
magnetic field to diameter (H — D) relation. This 
way, we can discuss SNR evolutionary properties by 
comparing theories on H with empirically extracted 
H — D relation. Statistical i.e. empirical study of H 
evolution nevertheless requires reliable data samples 
of SNRs in different types of interstellar medium. 

There is a number of ways to estimate H in SNRs. 
Unfortunately, few are reliable and even they are 
available only for a few well studied SNRs. The es- 
timates are made by measuring rotation measures 
or spectral line splitting. The estimates can also be 
extracted from radiation fluxes from different parts 
of the electromagnetic spectrum such as radio, X- 
rays or 7-rays. However, there is another problem in 
performing a statistical study of H based on these 
estimates. The data samples of SNRs are ballast 
by severe selection effects through the entire electro- 
magnetic spectrum. SNRs are mainly identified in 
the radio domain. Unlike optical, X-rays or 7-rays, 
radio waves are less influenced by absorption and 
scattering in the interstellar medium. Also, radio 
interferometers have the best resolution among all 
the other observational devices, which also helps in 
the detection of remnants. Large and reliable data 
samples are of crucial importance for a good and 
well-founded statistical study of empirical H — D 
relation. Today, this condition is partially fulfilled 
only by data samples in radio domain. The empirical 
studies of SNR properties are also severely influenced 
by the selection effects. It remains to be hoped that 
the observational instruments and techniques in the 
future will help us overcome this problem. 

The main purpose of this paper is to apply and 
discuss a method for the determination of H — D 
slope from radio luminosity at given frequency v to 
diameter [L v — D) correlation (£„ = L„/(D 2 ir 2 )), 
in SNR samples that show existence of such a corre- 
lation. The method is based on the energy equipar- 
tition assumption between magnetic field and rel- 
ativistic particles. It uses equations for revised 
equipartition calculation (REC). The equipartition 
calculation is the most commonly used manner of 
obtaining H estimates in valid radio SNR samples. 
The obtained H — D slope from the only reliable data 
sample of M82 SNRs is then compared with the slope 
arising from the theoretical models of SNR evolution. 
We can then argue whether or not M82 SNRs are in 
the equipartition state, and thereby give a contri- 
bution to the general evolutionary studies of SNRs. 
We also try to give estimates of the magnetic field 
strengths, particularly for SNRs in M82. In addition, 



we discuss the accuracy of magnetic field strength 
obtained under REC. This is done by comparing the 
values for H, obtained herein, with the more reli- 
able ones available in literature (found for few SNRs 
from Large Magellanic Cloud and our Galaxy). It is 
noteworthy that an SNR luminosity is mainly deter- 
mined by the density of environments in which SNR 
evolves. This is an important issue for the discussion 
of the influence of equipartition arguments on H. 

This paper is organized as follows: Section 2 
presents explanations of the required topics which 
are too broad to be mentioned later in the text. In 
Section 3 we describe and analyze the method and 
REC with its assumptions. Section 4 features a dis- 
cussion on the obtained results for H — D slope and 
magnetic field strength. There, we consider whether 
M82 SNRs are in the equipartition state or not. Fi- 
nally, the conclusions of this work are given in Sec- 
tion 5. 

2. THE H-D DEPENDENCE 

2.1. History of H — D Relation 

We assume that H — D relation can be written 
in the form: 

H(xD- 5 . (1) 

Historically, this form of the magnetic field evolution 
is used in all theoretical models that explain the syn- 
chrotron emission from SNRs. 

Shklovsky (1960) was the first to theoretically de- 
scribe the synchrotron emission from a spherically 
expanding nebula. He assumed that magnetic field 
structure remains unchanged during the expansion. 
Consequently, magnetic field flux is constant and 
H oc D~ 2 , where D is the diameter of the remnant. 
Lequeux (1962) applied Shklovsky's theory to model 
shell type remnants, which led to H oc (I x -D) -1 , 
where I oc D represents the thickness of the shell. 
Poveda & Woltjer (1968) and Kesteven (1968) also 
gave their contribution to the general model of shell 
type remnant. They assumed that H is gained with 
the compression of the interstellar medium magnetic 
field (leading to H = const) and that shell thick- 
ness remains constant during the expansion ( which 
leads to H oc D^ 1 ). Theoretical interpretation 
of SNR synchrotron emission by Duric & Seaquist 
(1986) used the magnetic field model with H oc D~ s , 
based on the work of Gull (1973) and Fedorenko 
(1983). According to the results of Gull, magnetic 
field is compressed and amplified in the convection 
zone, to finally gain enough strength to power the 
bright synchrotron emission. Fedorenko stated that 
1.5 < 8 < 2.0. Tagieva (2002) obtained H oc L>-° 8 , 
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by using the E oc £)~ 2 - 38 relation (Case & Bhat- 
tacharya 1998). However, this result should be taken 
with great reserve because the E — D relation from 
the work of Case & Bhattacharya is balast by severe 
selection effects (Urosevic et al. 2005). Also, Tagieva 
did not take into account the influence of the den- 
sity of environments in which the SNRs evolve. An 
interesting discussion about the magnetic field and 
the equipartition arguments for five Galactic SNRs, 
based on the results empirically obtained from an X- 
ray data, can be found in the work of Bamba et al. 
(2005). The predecessor of this paper is the work of 
Vukotic et al. (2006). 

2.2. Magnetic Field Calculation from Radio 

Synchrotron Luminosities 
The magnetic field is calculated from the follow- 
ing formula for synchrotron emission of relativistic 
electrons (Beck & Krause 2005, hereafter BK): 

L v = / tirfVc 2 (~f)n cfi ■ 

• E^^/2c 1 ) (1 -' )/2 H^ +1 ^ 2 . (2) 

We adjusted the formula from BK to suit our needs. 
Here, / is a fraction of the radio source volume 
occupied by the radiative shell. We assumed that 
/ = 0.25. This is consistent with SNRs having 
strong shocks where compression ratio is 4. How- 
ever, should this not be the case, variation of / will 
still not have any significant effect on values for H , 
because of the small value of exponent (7 + l)/2 in 
Eq. (2). Further, the total volume of SNR is desig- 
nated by V. Instead of spectral intensity along the 
radiation ray path (I u ) in BK, we used spectral lumi- 
nosity of the source, because the majority of sources 
in used data samples are seen almost as point-like 
sources, having only the flux density data integrated 
over the whole source available. According to BK 
this may lead to the overestimation of values for H. 
This effect is discussed further in Section 4. The rest 
is the same as in BK, 02(7) (in units erg -2 s _1 G _1 ) 
is identical to 05(7) in Pacholczyk (1970), n c , is 
the number density of cosmic ray electrons per unit 
energy interval for the normalization energy Eq, 
ci = 3e/(47rm3 c 5 ) = 6.26428 • 10 18 erg~ 2 s" 1 G"\ 
H± is the magnetic field strength in the sky plane, 
and finally 7 represents exponent in the cosmic ray 
power law energy spectrum (see Appendix A in BK) . 
Closer inspection of Eq. (2) shows that in order to 
calculate H from L v , some assumption regarding the 
relationship between H and n e ,o has to be made. 

2.3. Data Samples 
Currently, it seems that there is no better way to 
determine H by using only data on L v and spectral 



index a (7 = 2a + 1) than the equipartition or the 
minimum-energy assumption. This method is useful 
for SNR samples where all other data are lacking. 
However, Galactic SNR data samples are strongly 
biased by selection effects. The farther the object, 
the greater its brightness detection limit. The ex- 
tragalactic samples suffer from milder selection ef- 
fects. Their brightness detection limits (sensitivity 
lines) do not differ from one SNR to another be- 
cause all the SNRs in the sample are approximately 
at the same distance. In this study, we have re- 
lied on the only statistically reliable sample of SNRs 
from a nearby starburst galaxy M82 (Huang et al. 
1994). The equations that we used in calculating H 
are presented in Section 3. Inspection of those equa- 
tions shows that any H — D correlation requires the 
existence of L v — D correlation. If L v — D correla- 
tion does not exist, than it makes no sense to extract 
H — D relation from L v — D data. If SNR data sam- 
ples show no existence or poor L u — D correlation 
there are two possibilities: SNR luminosity does not 
evolve with the diameter, which is unlikely, or the 
sample is made of SNRs that evolve in different envi- 
ronments and is influenced by selection effects. This 
is explained in the next paragraph. 

In their work, Arbutina et al. (2004) showed 
that the best L v — D correlation exists for SNRs in 
M82. They also showed that some correlation exists 
for Galactic SNRs associated with large molecular 
clouds. Arbutina & Urosevic (2005) imply that the 
evolution of SNR radio surface brightness depends on 
the properties of the interstellar medium, primarily 
the density. They formed three SNR data samples 
from the existing ones (Galactic and extragalactic): 
the Galactic SNRs associated with large molecular 
clouds (GMC), oxygen-rich and Balmer-dominated 
SNRs. The main intent of Arbutina & Urosevic 
(2005) was to group SNRs by their properties, pri- 
marily the density of the interstellar medium in 
which they evolve (and also by SN type). They also 
argued that M82 sample is the best possible sam- 
ple that one can currently find. All SNRs from M82 
are likely to evolve in similar environment of dense 
molecular clouds. Consequently, they are very lumi- 
nous and being extragalactic, SNRs exhibit milder 
selection effects. The reliability of M82 sample is also 
discussed in Urosevic et al. (2005). By performing 
the Monte Carlo simulation, the authors showed that 
M82 sample is not severely affected by sensitivity se- 
lection effects, as in the case of other extragalactic 
samples (LMC, SMC, M31, M33). 

In this paper we have applied REC and calcu- 
lated 5 for M82 sample, as the best sample for statis- 
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tical study, and additionally we have analysed three 
samples from Arbutina & Urosevic (2005): GMC, 
oxygen-rich and Balmer dominated SNRs. We did 
not use the last three samples to calculate the slope 
5 because they are of poorer quality. Instead, we 
used them for checking the consistency of obtained 
H values with the global picture of SNRs evolution 
in different environments. Also, through literature 
search we found the magnetic field strengths for some 
SNRs from Table 2 and compared them with the 
values obtained in this paper. We searched the cat- 
alog of observational data on Galactic SNRs from 
Guseinov et al. (2003, 2004a, b) and papers avail- 
able on the Web-based Astrophysical Data Service 
( http:/ / adswww. harvard, edu/ ) . 3 

In calculation we have used the radio flux den- 
sity per unit frequency interval S v and radio spec- 
tral index a data (Table 2). These two properties 
are related as: 

S u =f3v~ a , (3) 

where /3 is the flux scale factor. The luminosity is 
calculated as L v = And 2 S v , where d is the distance 
to an SNR. In the case of extragalactic SNRs we 
assume that d is the same for all SNRs, and equal to 
the distance to the host galaxy. 

2.4. Magnetic Field and Relativistic Particles 

Since our studies are based on the radio syn- 
chrotron luminosity of SNRs, we can not treat 
magnetic field separately from relativistic particles. 
These two properties of an SNR are strongly coupled 
and it makes no sense to study them separately. 

As mentioned before, calculation of H from 
Equation (2) requires an assumption about n e ,o- 
This quantity also evolves with D. In Table 1 and 
Section 4.2 we present and discuss various assump- 
tions about n C fi(D) evolution and its effect on H(D) 
evolution ( assuming empirical L—D relation). Some 
of the n 0i o evolution patterns are only illustrative 
and are used for estimating the effect of different 
patterns on S. The pattern we used in our method 
to calculate H arises from the equipartition of en- 
ergies implying that energy densities stored in the 
magnetic field and relativistic particles are approxi- 
mately equal. The equipartition is widely used for H 
strength estimates, based purely on the radio data, 
in SNRs, galaxies, etc. It gives reasonably explain- 
able values for S and H . Taking all of this into ac- 
count we based our method on the equipartition of 
energies. 

3 ADS is NASA-funded project which maintains three bib- 
liographic databases containing more than 4.7 million records. 



Revised equipartition calculation (REC) used to 
calculate H is presented in detail in the work of BK. 
According to BK, REC gives better results than the 
classical equipartition calculation (CEC) presented 
by Pacholczyk (1970). 

2.5. Evolution of Magnetic Field in SNRs 

In this subsection we present the theoretical val- 
ues for 5 that characterize a particular SNR evolu- 
tion phase. These values, together with the ones ob- 
tained by our empirical method, are used in Section 

4 in the discussion of the most probable evolution 
scenarios for SNRs in M82. 

If SNRs are young, in early Sedov or free ex- 
pansion phase, they expand practically adiabatically, 
since radiative energy losses are negligible. Under 
the adiabatic expansion assumption i.e. conservation 
of energy in cosmic rays and magnetic field (-4 (W) = 
0), and equipartition conditions (u>cr = wh), where 
W is the total energy and the quantities u>cr and u>h 
are the energy densities of cosmic rays and magnetic 
field respectively, it follows that 6 = 1.5. Indeed: 

5W = |Woc^V)oc|(^), (4) 

4(W) =0=^ H otD- 3/2 , (5) 
dt 

where w is the total energy density. In conclusion, 
SNRs in the free expansion or early Sedov phase will 
have 6 = 1.5 if they are in the equipartition state. 

On the other hand, if SNRs are older, in the late 
Sedov or radiative phase, the value may be closer to 

5 = 1.25. The radiative phase is characterized by sig- 
nificant energy losses and SNR would later expand 
with the velocity v oc essure-driven snow- 
plow). If n Ci o oc ?ip,o oc n-^v (Bcrczhko & Volk 2004, 
hereafter BV), assuming equipartition H 2 oc n e ,o, S 
would be 5/4=1.25. The quantity n Pj o is the number 
density of cosmic ray protons per unit energy inter- 
val for the normalization energy Eq, and nn is the 
hydrogen number density. 

It is a general belief that, during the expansion, 
SNRs strongly amplify interstellar magnetic field. 
Two basic mechanisms of magnetic field amplifica- 
tion operate in SNRs. The first one is the Rayleigh- 
Taylor instability at the contact discontinuity be- 
tween the supernova ejecta and ISM swept by SNR 
forward shock. This scenario leads to 1.5 < 6 < 2 
(Fedorenko 1983) and is preferred in young SNRs. 
The second mechanism operates right behind the 
shock, where magnetic field is amplified by strongly 
excited magnetohydrodynamic waves. This is the 
probable mechanism for older remnants. 
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3. ANALYSIS AND RESULTS 

There are two most commonly used assumptions 
regarding the magnetic field and cosmic ray energy 
content: 1) the minimum of total energy stored in 
the particles and magnetic field, and 2) the equipar- 
tition between these energies. The minimum energy 
assumption gives 4/3 for the ratio of the energies 
stored in the particles and magnetic field, which 
is ~ 1. These two assumptions are thereby often 
treated as synonymous and both procedures are re- 
ferred to as the equipartition calculation. There are 
also two different methods for obtaining these two 
estimates: classical (Pacholczyk 1970) and revised 
(BK) equipartition, i.e. minimum-energy calcula- 
tion. We will only present the formulas that we have 
used in calculating H and reader is referred to the 
mentioned papers for a detailed treatment of the sub- 
ject. 

3.1. Classical Calculation 
Classical formulas are: 



4.5 2 / 7 (l + fc) 2/7 - 

. Cl2 2/r / -2/7 (£ , /2) -6/7 i 2/7 ) 

6 2 / 7 (l + fc) 2/7 - 

. Cl2 2/7/-2A(£>/2)- 6 / 7 L2/7. ( 7 ) 



In these expressions we have introduced the follow- 
ing quantities: k is the ratio of the energies of the 
heavy relativistic particles and relativistic electrons, 
C12 and C13 are functions which are weakly depen- 
dent on a and are tabulated by Pacholczyk (1970). 
The radio luminosity L integrated between radio 
synchrotron spectrum cutoff frequencies v\ and v 2 
is calculated as: 

,. ! y 2 =10 11 Hz 

L = Aird 2 / S v dv. (8) 

Jv 1 = W Hz 

Using Equation (3) we can eliminate (3 and obtain L. 
We used k — 40 which should be adequate for strong 
shocks in SNRs. Being luminous synchrotron emit- 
ters and having small linear diameters, SNRs from 
M82 arc likely to be young and have strong shocks, 
but their true nature is still a subject to debate. We 
obtained S from Equations (6) and (7) by replacing 
L with L v — D relation from Arbutina et al. (2004). 
Replacing L with L v does not have any noticeable 
effect on S. We also assumed that H depends on D 
only trough L or L v . Therefore, 



Ho,(D- 3 L u f 7 oc(D- 4A ) 



-4.4\2/7 



if L(D) oc D~ 1A (Arbutina et al. 2004). This gives 
5 = 1.26. To prove the assumptions in Equation (9), 
we have calculated L from equations (3) and (8), and 
H from Equation (7). Then we fit linear regression 
in log H- log D plane to obtain S = 1.26±0.08. This 
shows that 012(a) does not change with D, so it does 
not affect S, which is why we can calculate 5 directly 
from the slope s of the L v oc D~ s relation, 



(3 + *)-, 



(10) 



as in Eq. (9). Equations (6) and (7) differ to a 
constant, giving exactly the same 5. In the sequel we 
do not show results for minimum energy estimates of 
H. 

3.2. Revised Formulas 

The main revision of the classical formulas is in 
using K instead of k. The quantities K and k stand 
for the ratios of proton to electron number densities 
and energy densities, respectively. In the CEC, in- 
tegration of the radiation energy spectrum between 
fixed frequency limits is performed. As opposed to 
this in REC the integration is performed over the 
energy spectrum of relativistic particles. This gives 
more accurate results (see BK). 

The revised formulas are: 



f/cqp 

rev 



kxKA{>y,L v ,v,V,f,i) 

1 V(<*+3) 

C( 1 ,E 2 )(a + l) 



8TrKA( 7 ,L v ,v,V,f,i) 

l/(a+3) 



(11) 



C{l,E 2 ) 



(12) 



where 



^(7^2) = ^o 2 -{^(t) 



,7-2 



+ 



2-7 



(tr 2 -(t) 



7-2' 



1 , E 2 
- + In — 

2 E D 



for 7^2, (13) 
for 7 = 2, (14) 



and 



(9) 



A( 7 , L v , v, V, /, i) = J^l^l— (15) 
47rc 2 (7)£o 7 fVci{i) 

In the above equations the following quantities ap- 
pear: K is the ratio of proton-to-electron number 



G 



VUKOTIC ET AL. 



densities per particle energy interval for the normal- 
ization energy Eq, Ei presents the high-energy limit 
for the spectrum of cosmic ray particles. The spec- 
tral break at low energies for protons is designated 
as E p = 938.28 MeV = 1.5033 • lCT 3 crg and finally 
c±{i) is used to replace the projected field component 
H± with the total field H (see Appendix A in BK), 
with i being the projection angle. 

Equations (11) and (12) were originally taken 
from BK, with a few adjustments. To make equa- 
tions hold for 7 < 2 we used E 2 = 3 x 10 15 eV 
( Vink 2004) . Instead of K + 1 factor we used only 
K which is justified for proton-dominated plasma, 
and since the original formulas do not include the 
effect of possible synchrotron losses that affect the 
electron power law energy spectrum. Using K in- 
stead of K + 1 may provide an even better approxi- 
mation when taking into account synchrotron losses. 
To put it simple, it is as if there were almost no 
electrons in cosmic rays, and only protons remained. 
This can be justified by the fact that the protons 
are far more energetic than electrons and show less 
synchrotron losses. Such assumption does not have 
any significant effect on the values for H because of 
the l/(a+3) exponent in Equations (11) and (12). 
In this case, Equation (9) transforms into 

H oc (D- 3 L v ) 1/(a+3) oc ( D - 4 - 4 ) 1/(a+3) , (16) 
and Equation (10) becomes 

S=(3 + s)=L_. (17) 

In Eq. (16) we applied the L v — D correlation, to ob- 
tain S = 1.22, while fitting gives <5 = 1.19 ±0.08. For 
a we used an average spectral index of the whole 
sample (a = 0.6). The value for 5 from Equation 
(17), and the one obtained by fitting calculated val- 
ues for H using Equation (12), are almost identi- 
cal. The difference is negligible and we could have 
smoothly, as in CEC, calculate S from the slope of 
L v — D relation. 

In calculating H, we assumed that the mag- 
netic field in radiative shell of SNR is completely 
turbulent and has an isotropic angle distribution 
in three dimensions, giving c 4 = (2/3)' 7+1 ^ 4 (Ap- 
pendix A in BK). This is the best assumption to 
be made when the majority of SNRs are point-like 
sources, i.e. without maps for H. We also used 
K = {E p /E c ) h ' 1)/2 (Appendix A in BK), where 
E c = 511 keV = 8.187 10~ 7 erg designates the spec- 
tral break at low energies for electrons. The data 
for 21 SNRs from M82 from the work of Urosevic 



TABLE 1 
RESULTS FOR 5 



Direct 


Shklovsky (1960) (n e ,o oc £>-( 2 <*+ 3 )) 
Berezhko & Volk (2004) (n e ,o oc D^) 


0.125 
0.875 


Classical 


equipartition 




1.26 


Revised 


equipartition 




1.22 



et al. (2005), and the obtained values for H, are 
shown in Table 2. As it can be seen, the magnetic 
field strengths are up to 10 mG. Using L v /{A-KfV) 
in our formulas instead of I v /l (BK) could lead to 
an overestimation of the average field. Nevertheless, 
if the magnetic field is significantly overestimated it 
should not have a significant effect on the value for 
5. There is also a possibility that M82 remnants are 
pulsar driven wind nebulae (PWNe). Unlike shell 
type SNRs, PWNe have different mechanisms that 
maintain magnetic fields. Magnetic field strengths 
in PWNe are comparable with the ones we obtained 
from REC for M82 SNRs. This possibility is inves- 
tigated further in Section 3.4. 

3.3. Direct Derivation 

It is possible to derive <5 directly from Eq. (2) 
if there is an additional assumption concerning the 
evolution of n Ci0 with D. We consider models used 
by Shklovsky (1960), and the assumption of conser- 
vation of cosmic ray energy i.e. adiabatic expansion 
(e.g. BV). Respectively, these are 

n c , oc D-^+V (18) 

and 

n c , cx D-\ (19) 
Equation (2) together with L v — D relation gives 

/ n-4.4\ !/("+!) 

H oc . (20) 

V n c,0 / 

For an average spectral index a = 0.6 the results are 
presented in Table 1. Here, we found fitting unnec- 
essary because we already saw in Sections 3.1 and 
3.2 that the rest of the quantities from Equation 2 
do not change with D, at least not in a way to affects 
S. By using direct method we can only get values for 
H scaled to a constant because of proportionality of 
Equations (18) and (19). 
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3.4. Calculated and Literature-found H Values for 
GMC, Oxygen-rich and B aimer- dominated 
SNRs 

With a view to checking values obtained for H 
we performed the same REC on the SNRs associ- 
ated with large molecular clouds, oxygen-rich and 
Balmer dominated SNRs. According to Arbutina & 
Urosevic (2005), these SNRs form parallel tracks in 
radio surface brightness to diameter plane. If envi- 
ronmental density is higher we expect the SNR to 
be brighter. The implication is that SNRs with the 
same D should have different luminosities if environ- 
mental densities are different. According to Equa- 
tion (12), SNRs that evolve in a more dense environ- 
ment should also have stronger H than SNRs with 
the same diameter that evolve in a less dense envi- 
ronment. The data used and the obtained CEC and 
REC results for all groups of SNRs are presented in 
Table 2. 

Figure 1 presents a plot of all REC values from 
Table 2. It shows that SNRs in a more dense envi- 
ronment (M82, GMC, oxygen- rich) appear to form a 
track in H — D plane, while Balmer-dominated SNRs 
form another track that lies beneath the first one. 
Due to the dispersion and incompleteness of data 
samples, any statistical study of the tracks should be 
avoided, for now. We can, however, make some qual- 
itative conclusions. In Figure 1 we can see that REC 
does not change L v — D the evolution pattern. This 
is very convenient for the estimate of the reliability of 
H in M82 SNRs. From Figure 1 it is clear that H val- 
ues for SNRs in M82 seem consistent with the values 
for GMC and oxygen-rich remnants. They all evolve 
in a dense environment and accordingly may have a 
similar H — D evolution pattern. Their H values, ac- 
cording to Arbutina & Urosevic (2005), are different 
in comparison to the values for Balmer-dominated 
SNRs. This is because Balmer-dominated SNRs are 
likely to evolve in a low density environment. In the 
group that consists of Balmer-dominated, oxygen- 
rich and GMC SNRs, used in this work, we did not 
include PWNe, because REC is made for shell type 
SNRs. Accordingly, to avoid possible PWNe, we did 
not include SNRs with a < 0.4, which is the char- 
acteristic of PWNe (Gaensler & Slane 2006). From 
Figure 1 we can see that the most of SNRs in M82 
are, probably, not PWNe because they fit the evolu- 
tion pattern for SNRs in dense environments. In ad- 
dition, the higher spectral indices of the M82 SNRs 
(average a « 0.6; see Table 2) are not characteris- 
tic for PWNe. However, the possibility that at least 
some of these objects are PWNe should not be easily 
put aside. For now, we can only wait for the obser- 
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Fig. 1. The revised equipartition data in logH — logD 
plane. The SNRs are presented by: crosses (M82), open 
circles (oxygen-rich), open triangles (Galactic SNRs as- 
sociated with large molecular clouds), filled triangles 
(Balmer-dominated SNRs). 

vational instruments to advance, and for a possible 
detection of pulsars in M82. 

Table 2 also shows the best available literature- 
found values for H inferred from other methods, for 
Galactic and LMC SNRs. The agreement of these 
values with the values obtained from REC is another 
way to show the reliability of H estimates for SNRs 
in M82. This is one of the subjects discussed in Sec- 
tion 4. 

4. DISCUSSION 

4.1. Values Obtained for H 

Both, classical and revised equipartition calcula- 
tion contain various uncertainties and assumptions 
and as such, are of limited applicability (BK). Nev- 
ertheless, by performing CEC and REC we arrived 
at the conclusion that all of the imperfections do not 
have a noticeable effect on S, but could have a signifi- 
cant impact on the values for H. Inspection of Table 
2 shows that obtained H values are higher than those 
found in literature. Such overestimates are prob- 
ably due to replacement of I v with L v (BK). The 
assumptions regarding / and K in REC equations 
arc not of great importance because of the small 
l/(a + 3) exponent. Due to the L v — ► I v replace- 
ment, the amount of overestimate is strongly affected 
by SNR morphology and consistently shows consid- 
erable variations from one SNR to another (Table 
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2). The morphology variations should not depend 
on diameter, which means that overestimates of H 
are mainly arising from morphology related factors, 
and they should only make data scattering in H — D 
plane without affecting 5. Table 2 also shows that an 
average overestimate by a factor of 2 can be adopted. 
Coupled with explanation of Figure 1 (Section 3.4), 
this shows that H values for SNRs in M82 are esti- 
mated reliably to an order of magnitude. This means 
that M82 does contain SNRs with magnetic fields of 
up to 10~ 2 G. However, this should be taken with 
reserve because of the possibility that some SNRs in 
M82 are perhaps PWNe. 

REC used in this paper thus gives reliable esti- 
mates accurate to an order of magnitude. This is of 
small significance in studies of nearby, well resolved 
SNRs with data from all parts of electromagnetic 
spectrum, but may be of great applicability in statis- 
tical and empirical studies of SNRs residing in other 
galaxies, that are unresolved and often have only ra- 
dio data available. As already mentioned, Galac- 
tic SNR samples are strongly influenced by selection 
effects and can not be used in the statistical and 
empirical studies of SNRs evolution properties. For 
now, the only SNR samples that can be used for reli- 
able statistical and empirical studies reside in other 
galaxies. With these samples, the obtained values for 
H will be probably overestimated by a factor of 2, 
but accurate to an order of magnitude, as in this pa- 
per. In the next section we discuss the results on the 
magnetic field evolution obtained when our method 
is applied to SNRs in M82. This should illustrate 
how the method can be used for getting closer in- 
sight into the SNRs evolution properties, i.e. SNRs 
evolution phases, and how it can be used to check 
the validity of the equipartition assumption. 

4.2. Magnetic Field Evolution of SNRs in M82 

If the sample is statistically reliable, the obtained 
H may be overestimated, but chosen REC param- 
eter values should not have a significant effect on 
6. The difference between S obtained from classi- 
cal and revised methods is mainly due to the ex- 
ponents in equations (7) and (12). These expo- 
nents will be equal for an average spectral index 
a = 0.5. For SNRs in M82 a = 0.6 is used, 
and therefore we obtain slightly different slopes in 
H — D plane. In the work of Bcrkhuijscn (1986), 
the author implies that a could depend on the den- 
sity of the ISM in which the SNRs evolve as: a = 
(0.075 ± 0.024) log n + (0.538 ± 0.012), where n is 
the density of the ISM. This means, according to 
Eq. (17), that the lower track in Figure 1, that con- 



sist of SNRs evoloving in the small density environ- 
ments, should have somewhat shallower slope when 
compared with the track above (large density envi- 
ronments). However, cosidering Eq. (17) and just 
mentioned relation, it is clear that for typical val- 
ues of a and no, there will be no significant effect 
on S. Consequently, the tracks from Fig. 1 should 
be considered as parallel. Taking all of above into 
account, we conclude that 5 is strongly affected by 
the assumptions regarding n c .o- "Directly" obtained 
values for 5 of 0.125 and 0.875 (Table 1) arc only 
illustrative. Shklovsky's model have a rather histor- 
ical meaning, since no additional particle accelera- 
tion (by the shock) during evolution is assumed ( 
besides the initial acceleration in supernova explo- 
sion). This leaves us with the equipartition as our 
best assumption. 

Table 1 shows that the equipartition arguments 
combined with the possible L u — D dependence give 

5 « 1.2. This value is slightly lower than theoretical 
value <5=1.5 obtained under equipartition and adia- 
batic approximations (Section 2.5). If SNRs in M82 
are young, in early Sedov or free expansion phase, 
this difference can be explained by the sensitivity se- 
lection effects related to the M82 sample. The Monte 
Carlo simulations in Urosevic et al. (2005) show that 
the measured slopes of extragalactic surface bright- 
ness to diameter (Ej,— D) relations are shallower due 
to the sensitivity selection effects. Therefore, the ap- 
parent £„ — D ( and L v — D) slope for M82 is lower 
than the real slope. The lower L v — D slope gives 
lower S. This means that equipartition arguments 
for the SNRs in M82 sample may still be applicable, 
whereas a small difference between the theoretical 
and empirical 6 can be ascribed to selection effects. 

On the other hand, S = 1.2 might indicate that 
not all SNRs from M82 sample are in the equiparti- 
tion stadium. If, for example, the larger ones are in 
the late Sedov phase where magnetic field remains 
constant (BV), empirical S would be a compromise 
between values and 1.5. The evolutionary status 
of SNRs remains a great uncertainty. The SNRs in 
M82 may be in the free expansion, as well as in the 
Sedov's, or even in the radiative phase. Chevalier 

6 Fransson (2001) proposed that M82 SNRs may 
be in the radiative phase because they evolve in a 
very dense environment. In this case, S may be 1.25, 
close to the empirical value. As the previous ones, 
this scenario too, remains uncertain. 

4.3. Interstellar Magnetic Field in M82 
Condon (1992) estimated field strength in M82 to 
be H w 100 /iG from classical minimum energy cal- 
culation, considering that central emitting region of 
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M82 is 30" x 10" and probably 0.5 kpc thick. Har- 
grave (1974) estimated the central emitting region 
in M82 to be 50" x 15". Using revised equipartition 
we estimated the value of « 190 [iG for the average 
interstellar magnetic field in the central emitting re- 
gion of M82 using the data Si. 4 gh z = 8.2 Jy and 
a = 0.68 from Klein et al. (1988). We assumed 
that / = 1 and that M82 radiates mainly from its 
central region of ~ 500 pc in diameter. This esti- 
mate is rough and should be taken with some re- 
serve. Such ISM magnetic field strength is among 
the highest field strengths when compared to other 
galaxies. This, however, may imply that M82 cen- 
tral region contains interstellar matter made of very 
dense molecular clouds. This is consistent with the 
high values of H in M82 SNRs, supporting the pos- 
sibility that their luminous synchrotron emission is 
mainly due to very dense environments and not due 
to pulsar driven wind nebulae. 

The values of up to 10 mG for H in M82 SNRs 
however imply that the magnetic field is strongly am- 
plified from the average ISM values of 100 — 200 ^G. 

5. CONCLUSIONS 
In this paper we presented and discussed a 
method for the determination of the magnetic field 
evolution pattern in SNRs only from the radio lu- 
minosity data samples. Such samples are the only 
available for the statistical and empirical studies of 
SNR evolution properties. The best sample, for now, 
consists of SNRs in M82, since these remnants seem 
to evolve in similar environment and share similar 
properties, and are not severely influenced by selec- 
tion effects. 

In order to calculate H from REC we were forced 
to make some assumptions. The only significant 
effect on values for H, regarding the assumptions, 
comes from replacing I v in REC formulas from BK 
with L„, which is done in order to apply REC on 
practically point-like sources. The other assump- 
tions are less important because of the small expo- 
nent in REC equations. Obtained under equiparti- 
tion assumption, 8 is a direct consequence of L v — D 
slope and has reasonable theoretical explanation. All 
assumptions do not change the evolutionary picture 
from L v — D plane. This means that our empiri- 
cal estimate of 5 is likely to be reliably determined. 
When compared with the more reliable values found 
in literature, the obtained H values appear to be 
overestimated approximately by a factor of 2. We 
conclude that H values for all SNRs, even the ones 
from M82, are accurate to an order of magnitude. 

To answer whether or not M82 SNRs are in 
equipartition state we have compared empirical 8 ob- 



tained by our method with the theoretical values. 
The empirically obtained <5 from L v — D correlation 
under the equipartition assumption is probably the- 
oretically explainable by the following two scenarios: 

(i) The slight difference between the theoretically 
derived H—D slope (8 = 1.5) under the adiabatic ap- 
proximation and the equipartition assumption, and 
the slope obtained in this paper using the empiri- 
cal L v — D correlation and REC (8 « 1.2) can be 
explained by the sensitivity selection effects which 
affected the sample of SNRs in M82. In this way, 
the starting assumption concerning the approxima- 
tive equipartition between the energy stored in the 
relativistic particles and in the magnetic field, could 
be justified. Therefore, we can conclude that SNRs 
in the M82 sample are probably close to the equipar- 
tition state. 

(ii) Finally, equipartition conditions may not be 
fulfilled for all remnants. If, for instance, they are 
in different stages of evolution, 8 may be between 
and 1.5. 

If SNRs are in the adiabatic phase, the most 
probable explanation for the lower empirically ob- 
tained value for 8 are the sensitivity selection effects 
in the M82 sample, perhaps in combination with 
slight deviation from equipartition, but the problem 
is the unresolved evolutionary status of M82 SNRs. 
Additional observations of SNRs in nearby starburst 
galaxies are needed for any firmer conclusions to be 
made. 
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TABLE 2 
SNRS DATA a AND RESULTS 



Catalog 


Other 


Type 1 


D 


Si 


a 


Distance 




rrcqp 
1 1 rev 




name 


name 






flux density 
at 1 GHz 






class, 
equip. 


rev. 
equip. 


literature 








(pc) 


(mJy) 




(kpc) 


(G) 


(G) 


(G) 


M82 39.1+57.4 




MC 


0.9 


8.28 


0.50 


3.9 x 10 3 


6.03E-03 


8.76E-03 




M82 39.4+56.1 




MC 


3.23 


4.25 


0.58 


3.9 x 10 3 


1.68E-03 


2.10E-03 




M82 39.6+53.4 




MC 


2.65 


2.68 


0.45 


3.9 x 10 3 


1.74E-03 


2.96E-03 




M82 40.6+56.1 




MC 


3.02 


4.97 


0.72 


3.9 x 10 3 


1.94E-03 


2.24E-03 




M82 40.7+55.1 




MC 


1.93 


15.56 


0.58 


3.9 x 10 3 


3.78E-03 


4.64E-03 




M82 41.3+59.6 




MC 


1.02 


6.19 


0.52 


3.9 x 10 3 


4.99E-03 


6.85E-03 




M82 42.7+55.7 




MC 


4.30 


6.10 


0.71 


3.9 x 10 3 


1.51E-03 


1.78E-03 




M82 42.8+61.3 




MC 


1.97 


3.58 


0.63 


3.9 x 10 3 


2.47E-03 


2.92E-03 




M82 43.2+58.4 




MC 


1.05 


12.61 


0.66 


3.9 x 10 3 


6.11E-03 


6.83E-03 




M82 43.3+59.2 




MC 


0.60 


29.54 


0.68 


3.9 x 10 3 


1.27E-02 


1.35E-02 




M82 44.3+59.3 




MC 


1.96 


5.46 


0.64 


3.9 x 10 3 


2.80E-03 


3.27E-03 




M82 44.5+58.2 




MC 


2.25 


3.55 


0.50 


3.9 x 10 3 


2.16E-03 


3.13E-03 




M82 45.2+61.3 




MC 


1.12 


19.54 


0.67 


3.9 x 10 3 


6.58E-03 


7.28E-03 




M82 45.3+65.2 




MC 


2.05 


5.80 


0.82 


3.9 x 10 3 


2.96E-03 


3.32E-03 




M82 45.4+67.4 




MC 


2.23 


5.01 


0.67 


3.9 x 10 3 


2.47E-03 


2.86E-03 




M82 45.8+65.3 




MC 


2.13 


3.74 


0.46 


3.9 x 10 3 


2.30E-03 


3.79E-03 




M82 45.9+63.9 




MC 


2.22 


4.25 


0.41 


3.9 x 10 3 


2.32E-03 


4.70E-03 




M82 46.5+63.9 




MC 


1.39 


6.93 


0.74 


3.9 x 10 3 


4.18E-03 


4.60E-03 




M82 46.7+67.0 




MC 


2.95 


4.39 


0.76 


3.9 x 10 3 


1.94E-03 


2.25E-03 




M82 41.9+58.0 




MC 


0.52 


154.96 


0.75 


3.9 x 10 3 


2.38E-02 


2.32E-02 




M82 44.0+59.6 




MC 


0.79 


54.89 


0.48 


3.9 x 10 3 


1.16E-02 


1.80E-02 




G 111.7-2.1 


Cas A 





4.9 


2720 x 10 3 


0.77 


3.4 


1.02E-03 


1.23E-03 


5.5E-04 b 


G 260.4-3.4 


Pup A 





35.2 


130 x 10 3 


0.5 


2.2 


5.73E-05 


8.32E-05 




LMC 0525-69.6 


N132 D 


o 


25 


5800 


0.7 


55 


2.07E-04 


2.71E-04 


< 4E-05 C 


SMC 0103-72.6 







55 


250 


0.5 


65 


4.53E-05 


6.58E-05 




NGC 4449 







0.6 


20 


0.75 


4200 


1.22E-02 


1.25E-02 




G 42.8+0.6 




MC 


76.8 


3 x 10 3 


0.5 


11 


2.51E-05 


3.64E-05 




G 78.2+2.1 


7 Cygni 


MC 


20.9 


340 x 10 3 


0.5 


1.2 


8.34E-05 


1.21E-04 




G 84.2-0.8 




MC 


23.6 


11 x 10 3 


0.5 


4.5 


6.00E-05 


8.71E-05 




G 89.0+4.7 


HB 21 


MC 


24.2 


220 x 10 3 


0.4 


0.8 


5.20E-05 


9.87E-05 




G 132.7+1.3 


HB 3 


MC 


51.2 


45 x 10 3 


0.6 


2.2 


3.10E-05 


4.24E-05 




G 166.2+2.5 


OA 184 


MC 


183.8 


11 x 10 3 


0.57 


8 


1.43E-05 


2.08E-05 




G 309.8+0.0 




MC 


23 


17 x 10 3 


0.5 


3.6 


6.11E-05 


8.88E-05 




G 315.4-2.3 


MSH 14-63 


MC 


28.1 


49 x 10 3 


0.6 


2.3 


5.45E-05 


7.34E-05 




G 349.7+0.2 




MC 


8.7 


20 x 10 3 


0.5 


14.8 


3.3E-04 


4.80E-04 


3.5E-04 d 


G 4.5+6.8 


Kepler 


B 


2.4 


19 x 10 3 


0.64 


2.9 


3.95E-04 


4.97E-04 


2.15E-04 b 


G 120.1+1.4 


Tycho 


B 


5 


56 x 10 3 


0.61 


2.3 


2.49E-04 


3.21E-04 


3E-04 b 


G 327.6+14.6 


SN 1006 


B 


19 


19 x 10 3 


0.6 


2.2 


5.67E-05 


7.62E-05 


1.6E-04 b 


LMC 0505-67.9 


DEM L71 


B 


19 


9 


0.5 


55 


3.96E-05 


5.75E-05 




LMC 0509-68.7 


N103 B 


B 


7 


1100 


0.6 


55 


3.72E-04 


4.75E-04 




LMC 0509-67.5 




B 


7 


70 


0.5 


55 


1.68E-04 


2.43E-04 




LMC 0519-69.0 




B 


8 


150 


0.5 


55 


1.86E-04 


2.70E-04 




LMC 0548-70.4 




B 


25 


100 


0.6 


55 


6.29E-05 


8.44E-05 




SMC 0104-72.3 




B 


29 


12 


0.5 


65 


3.29E-05 


4.78E-05 





Notes: a M82 data are taken from Table A.l in Uroscvic ct al. (2005) with Si being scaled from 1.4 to 1 GHz. The rest of the 
used data are same as in papers of Arbutina et al. (2004) and Arbutina & Uroscvic (2005), with data of Galactic MC SNRs being 
updated for distances from Green (2004); b V61k et al. (2005); c Dickcl & Milne (1995); d Brogan et al. (2000); 1 MC - Associated 
with giant molecular clouds, O - Oxygen-rich, B - Balmcr-dominatcd. 
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